Magnetic screening in proximity effect Josephson-junction arrays 
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Abstract. - The modulation with magnetic field of the sheet inductance measured on proximity 
effect Josephson-junction arrays (JJAs) is progressively vanishing on lowering the temperature, 
leading to a low temperature field-independent response. This behaviour is consistent with the 
decrease of the two-dimensional penetration length below the lattice parameter. Low temperature 
data are quantitatively compared with theoretical predictions based on the XY model in absence 
of thermal fiuctuations. The results show that the description of a JJA within the XY model is 
incomplete and the system is put well beyond the weak screening limit which is usually assumed 
in order to invoke the well known frustrated XY model describing classical Josephson-junction 
arrays. 



Introduction. — Proximity effect Josephson-junction 
arrays (JJAs) exposed to a perpendicular magnetic field B 
are well known to be, under some conditions discussed be- 
low, a good physical realization of the frustrated XY model 
[1] and so have been studied for many years as model sys- 
tems, and still reveal some interesting physics [2]. The 
description of JJAs in the framework of the frustrated XY 
model remains valid in a broad temperature range below 
the superconducting-to-normal transition temperature Te- 
la the case of proximity effect SNS JJAs, quantum fluc- 
tuations are irrelevant even at very low temperature. In- 
deed, the Josephson coupling energy Ej = (PqIc/^-k (0o 
the flux quantum and Ic the single junction critical cur- 
, rent) is much larger than the Coulomb charging energy 
' Ec = e*^/2C (e* = 2e the free Cooper pair charge and 
C the junction capacitance), and the normal state resis- 
tance is much smaller than the quantum resistance [3] 
(roughly 7 orders of magnitude smaller in our systems). 
These two criteria ensure that our JJAs remain in the 
classical regime for all T < Tc- However, when lowering 
the temperature magnetic interactions (between screening 
currents of neighbouring cells) put the system beyond the 
pure XY model. As a matter of fact, the system expe- 
riences a gradual crossover from a high temperature XY 
weak screening regime where vortex ground states satisfy 
the fluxoid quantisation in all plaquettes [1], and local 
magnetic fleld equals the applied field B, to a low temper- 



ature strong screening regime with quantised fiux where B 
penetrates the JJA in a similar way to how it penetrates a 
superconducting wire network creating an Abrikosov-like 
lattice. In the weak screening limit, the frustration pa- 
rameter / is defined by the magnetic flux (j) threading an 
elementary cell in units of the flux quantum (j)^] f — (j)/4>o. 
While lowering the temperature, the JJA is gradually put 
beyond the weak screening limit when supercurrents flow- 
ing around each lattice cell are no longer negligible. At 
very low temperatures, when an elementary cell of the 
array is able to fully screen a flux quantum, any connec- 
tion to the XY model is lost since the JJA behaves as a 
multiply-connected superconductor. 

Although their existence is well known, the effects of 
screening currents on the behaviour of JJAs have been 
investigated mainly theoretically, in the weak screening 
regime (as a correction of the frustrated XY model) [4] 
and beyond [5,6]. This paper focuses on the low temper- 
ature regime of JJAs on a periodic lattice with hexagonal 
symmetry called a dice lattice (see Fig. 1 and description 
below) where the measured sheet inductance is indepen- 
dent of B. In this regime, three different phenomena push 
our JJA beyond the XY model. Firstly, the effective pene- 
tration length becomes shorter than the lattice parameter 
and the magnetic energy associated with screening cur- 
rents dominates over the Josephson coupling energy as 
far as the magnetic field response is concerned. Secondly, 
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thermal fluctuations gradually vanish and vortex mobility 
is accordingly reduced. Finally, on a smaller size scale, 
the Josephson penetration length becomes smaller than 
the size of a Josephson junction which results in a non- 
unique phase difference across the junction width. All this 
leads to a low temperature regime dominated by magnetic 
screening effects which efficiently suppress the field mod- 
ulation of the measured sheet inductance. This scenario 
is confirmed by comparing low temperature sheet induc- 
tance data with the same quantity calculated within the 
frustrated XY model in absence of thermal fluctuations, 
i.e. in the ground states. Thus, when the temperature 
is lowered JJAs can easily reach a regime where the XY 
model is not sufficient to describe their behaviour. 
In addition, for the first time data taken using the two- 
coil mutual inductance technique [7], as well as the re- 
lated data treatment, are confirmed by comparing Ic ex- 
tracted from inductive measurements with the same quan- 
tity measured by a conventional four probe measurement 
technique. 

Measurement technique. — Arrays of proximity ef- 
fect Pb-Cu-Pb junctions on a dice lattice (see upper inset 
to fig. [U were probed using a SQUID-operated two-coil 
mutual inductance technique [7] that allows for measure- 
ments of the screening properties of superconducting wire 
networks [8] , JJAs [2] , as well as high temperature super- 
conducting thin films [9] . Here the macroscopic measured 
quantity is the inverse sheet inductance (or inverse 
inductance "per square" [10]) that is inferred from the 
array's linear sheet conductance [7], and measures the 
degree of superconducting phase coherence in the sample. 
The inductive technique has a sensitivity threshold of the 
order of lOpH. 

For comparisons with theoretical predictions based on 
the XY model and with samples with different coupling 
constants, i.e. different geometrical parameters, the tem- 
perature is scaled according to the thermal energy fc^T 
and Ej\ the so-called reduced temperature is r = kBT/Ej. 
At zero field and at temperatures well below Tc the sheet 
inductance of a regular array L cx Lj [10] the single junc- 
tion inductance, and the numerical factor depends only on 
the lattice geometry. For the dice lattice [11] L^^{T) — 
{2/y/i)L-\T) and L-\T) = (2/V3)(2V</>o)/c(T) [12]. 
Close to Tc, the critical current Ic is extrapolated from 
low temperature data [13]. 

For the first time inductive measurements performed us- 
ing the two-coil mutual inductance technique [7] , and the 
successive data treatment, are independently and directly 
confirmed by comparing the temperature dependence of 
the single junction critical current Ic (at zero field) ex- 
tracted from inductive measurements on a JJA with the 
same quantity measured by four probe measurements fol- 
lowing the procedure detailed in ref. [14] where Ic is de- 
fined as the current that produces a maximum in the dy- 
namic resistance {dV/dl vs /). Since the small normal 
state sheet resistance of proximity effect JJAs (of the or- 
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Fig. 1: (colour on-line) Low temperature dependence of the 
single junction critical current Ic extracted from resistive mea- 
surements (red plain symbols) and measured inductively (black 
empty symbols) at zero field. Upper inset: SEM picture of a 
portion of a 1 x Icm^ JJA on dice lattice. The array consists 
of ~ 10^ star-shaped (6-fold and 3-fold) superconducting lead 
islands coupled by proximity effect junctions in the underly- 
ing normal copper layer. The junctions have nominal width 
w = 2fim and length I = 1/im. An elementary rhombic cell 
illustrated with black line has a side a = 8fim. Lower inset: 
photograph taken with an optical microscope of modified sam- 
ple for resistive dV/dl measurements. The current / is injected 
on copper contacts at bottom into the central zigzag strip, and 
voltage V -I- /V — is measured on copper contacts at the top 
ends of the strip. The strip consists of 18 branches of width 
about containing 48 parallel junctions. The branches 

are equally spaced by 100/im and go back-and-forth across an 
area slightly smaller than the initial JJA. The strip is about 
12cm long. 

der of a mSl) and the high current that one would have 
to apply to approx. 1000 parallel channels of junctions in 
the Icm^ sample, and leading to a noticeable increase of 
temperature in the array, the sample geometry was mod- 
ified after all inductive measurements were done. A new 
pattern is obtained by an additional optical lithography 
process with the photoresist layer spread out on the top 
of the pre-existing JJA, and part of the superconduct- 
ing lead islands were removed using Ar-ion milling. The 
resulting pattern was a thin (48 parallel junctions) and 
12cm long back-and-forth strip across the initial array, as 
shown in the lower inset of fig. [1] The pre-existing copper 
layer all around the JJA is left on the substrate and used 
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for current and voltage contacts, as well as to dissipate 
the heat due to the measurement current (not shown in 
fig- El) ■ This new pattern allows for a lower applied current 
and the measured voltage drop across the strip is much 
greater than it would be across the initial array. Particu- 
lar care during the modification process was taken to avoid 
residual lead islands in-between the branches of the strip 
that could create electric shortcuts between the branches. 
Fig. [1] shows low temperature values of Ic obtained from 
both inductive measurements on the original JJA and re- 
sistive measurements on the modified sample. The error 
in /c obtained by resistive measurements is due to the res- 
olution in measuring the current and the related position 
of the maximum in dV/dl. The measurements are limited 
to temperatures below 4.2K because the sample was im- 
mersed in the helium bath of the cryostat to avoid local 
heating of the junctions which could be induced by the 
applied current. Insight into the effects of the reduction 
in the size of the array can be provided by a comparison 
of our experimental conditions with previous results ob- 
tained by numerical simulations of the resistively shunted 
Josephson-junction model at zero temperature (T = 0) 
[5], thus neglecting finite temperature effects. At zero field 
and at T = the dynamics of a JJA reduces to that of 
a set of uncoupled channels of single junctions along the 
current direction [5]. In our experiments the ambient 
magnetic field is reduced by a combination of mu-mctal 
and superconducting screens and at nominally zero field 
the residual frustration is / = 10"'^ {B = 0.4mG), which 
we can reasonably assume small enough to compare our 
data with results at / = [5]. A frustration parameter 
f = 1, i.e. one quantum of flux per cell, is achieved with 
a perpendicular magnetic field B = 360mG. When edge 
fields induced by the applied external current are taken 
into account and in the case of a 32x32 array [4] (com- 
pared to our modified sample with 48 parallel junctions), 
Ic is reduced by less than 10% when the field penetration 
length is of the order of the lattice parameter, i.e. in the 
strong screening regime (see discussion below). In the fol- 
lowing we show that in our JJAs such a strong screening 
regime is achieved for temperatures T < 3K. Thus we do 
not expect our resistive measurements to be significantly 
affected by the reduced size of the sample, except possibly 
for the few measurements below 3K in fig. [T] The data 
presented in Fig. [1] show the validity of the measurement 
technique, as well as the procedure to extract the critical 
current. 

Results and discussion. — In the weak screen- 
ing regime, due to the nature of the coupling in JJAs 
[1], the frustration / can induce a pronounced mod- 
ulation of the single junction inductance Lj{T) = 
{4>o/2tt) / {Ic{T)cos{9)) [15], hence the measured sheet in- 
ductance L{T), where 9 is the gauge invariant phase dif- 
ference related to / by the fluxoid quantisation [1,2]. 
Fig. [2] shows the response of a JJA while sweeping B in 
the reduced interval < / < 1/2 [1]. 
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Fig. 2: (colour on-line) Inverse magnetoinductance L~^{f) 
isotherms measured in the temperature range between r = 0.22 
(T=5K) and r = 0.004 (T=3K) with an excitation frequency 
u;/27r=16Hz. 

Well defined structures are measured at low excitation 
frequency showing peaks in L~^(f) which reveal enhanced 
superconducting phase coherence for some rational values 
of frustration (/ = 0,1/6,2/9,1/3). The / = 1/3 state 
has been the subject of previous work which shows that 
a strongly enhanced superconducting phase coherence, as 
for / = 1/3, does not necessarily imply that the vortex 
pattern is ordered [16, 17], even if ordering in the vortex 
pattern is normally associated with the commensurability 
between the vortex lattice induced by B and the array 
geometry, and resulting in a strong phase coherence as 
shown in fig. [2] for / = 1/3. The fully frustrated state 
(/ = 1/2), which on the contrary shows a pronounced 
depressed response in the intermediate temperature 
range, was also investigated [18, 19] and for both / = 1/3 
and / = 1/2 states the magnetic interaction between 
screening currents circulating in neighbouring cells of 
the array has been considered, but still assuming a 
uniform frustration over the array. Starting from the 
hottest isotherms and lowering the temperature, L^^{f) 
first shows the appearance of well defined structures at 
/ = 0,1/6,1/3 which grow up until some temperature 
T w 0.1. While further lowering the temperature the 
structures progressively disappear as shown by the coldest 
isotherms in fig. [5J 

Clues about the very low temperature response (r < 
10^^) come from the temperature behaviour of the array 
penetration length. For a thin superconducting film of 
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Fig. 3: (colour on-line) Low temperature dependence of the 
effective penetration length A in unit of the lattice parameter 
a = 8/im (see upper inset to Fig.[T| for three frustration values 
/ = 0,1/3,1/2. Inset: energy ratio E^/Ej as a function of 
reduced temperature r. 



thickness d and exposed to B the screening distance is 
given by a thickness-dependent penetration length [20], 
the bulk material (3D) penetration depth A being rescaled 
according to the thickness d, the effective 2D penetra- 
tion length is A = 2X^/d. In the case of an array, A 
can be expressed in terms of Ic [12], or as mentioned 
above as a function of the measured sheet inductance; 
A = 2{2/\/S)L/ fiQ for a dice lattice. Fig. [3] shows the 
low temperature dependence of A expressed in unit of the 
lattice parameter a for three values of frustration: / = 0, 
/ = 1/3, and / = 1/2. 

At very low temperature, below r « 3 • 10"'^, the 
field dependence of A (i.e. the superconducting phase 
coherence) is suppressed in good agreement with the 
magnetoinductance isotherms of fig. [H Most interestingly, 
this temperature is very close to the temperature at which 
A — a, the lattice parameter. In other words, for A close 
to or smaller than the lattice parameter the sample no 
longer behaves as a JJA with characteristic size a. This 
observation suggests the presence of strong magnetic 
screening. More information is obtained by comparing the 
magnetic energy Em stored in each supercurrent Ig loop 
around the lattice plaquettes with the Josephson energy 
Ej. For the sake of simplicity E^ is overestimated by 
considering the maximum supercurrent, i.e. the critical 
current Ic; Em = ^^Is ^ ^^^c = ^m- The geometrical 
inductance C = 26pH of a rhombic loop (see upper 
inset to fig. [T]) is calculated using the numerical results 
of [21] for a wire with rectangular cross section. The 
ratio E*^/Ej shown in the inset of fig. [3] shows that the 
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Fig. 4: Josephson penetration length Aj as a function of re- 
duced temperature r. The dashed line indicates the size of a 
juntion (width w — 2/im). 

magnetic energy is dominating over the Josephson energy 
for temperatures r < 5 • 10"'^. 

Screening effects also appear locally on each Joseph- 
son junction, inside which Meissner screening fields 
are responsible for inhomogeneous distribution of the 
Josephson currents. Quantitatively this phenomenon is 
negligible for small junctions, i.e. when the junction 
width w ^ Xj the Josephson penetration length given by 
A,/ = y^h/ {2e^QlJc) [22], where Jc is the cross section 
critical current density and I is the junction length 
(see fig. [TJ. Fig. 2] shows the temperature dependence 
of Aj which is smaller than the junction width w for 
T < 4 • lO^'^. Thus, at lower temperature the current dis- 
tribution inside the junctions is no longer homogeneous. 
As a consequence, the phase difference across the junction 
width is no longer unique and the JJA cannot be mapped 
onto the XY model. 

There exist two different mechanisms for the suppres- 
sion of the magnetic field dependence of the sheet induc- 
tance measured at finite frequency while the temperature 
is lowered: magnetic screening induced by superconduct- 
ing currents in the array and the vanishing of thermal 
fluctuations leading to the freezing of the vortex diffusion 
between array cells. Evidence for the predominance of the 
flrst mechanism at the lowest temperatures comes from 
a comparison of the measured sheet inductance with the 
same quantity calculated in the XY regime and in ab- 
sence of thermal fluctuations, i.e. the ground state value 
of L^^{f) calculated within the frustrated XY model at 
selected frustration states / = 0,1/3,1/2. At very low 
temperature where thermal fluctuations can be neglected. 
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the JJA is equivalent to a network; of inductances [23] 
Lij {qij ) and the phase variables qij are known in the frus- 
trated ground states at / = 1/3 [17] and at / = 1/2 [19]. 
Each ground state is characterised by a sheet inductance 
L{f) that can be calculated from the frustrated inductance 
network L^j which is first transformed into an anisotropic 
triangular lattice L'^j by the triangle-star transformation 
[15]. L{f) is then calculated applying Kirchhoff's laws to 
the elementary cell of the transformed network L'^^ . 

Unfrustrated state / = 0. All bonds of the trans- 
formed triangular lattice have the same inductance value, 
hence the isotropic sheet inductance of the dice lattice 
L(0) = {■\/3/2)Lj with Lj the single junction inductance 
at zero field. 

Frustrated state f — 1/3. Here we consider the 
(honeycomb) ground state which is dominating on a 
short scale over the possible vortex configurations [17]. 
The transformed sheet inductance is isotropic; L(l/3) = 
(4\/3/5)L.7 = 1.39Lj. 

Fully frustrated state f = 1/2. Applying the triangle- 
star transformation to the four periodic ground states 
[19], we obtain two different configurations each with two 
inductance components. Thus, in the ground state of 
the fully frustrated XY model the JJA behaves as a two 
component inductance network (Lx,Ly). Moreover, due 
to the inefficiency of the degeneracy removal mechanism 
and the prominence of the finite size effects [19] the 
vortex pattern can be assumed as disordered. In a 2D 
disordered system with two phases equally distributed, 
the conductivity is given by the geometric mean of both 
phase conductivities [24]. Although in principle one 
should then consider the geometric mean of {Lx,Ly), 
the discussion on the choice of the mean value, i.e. the 
arithmetic or the geometric mean, is irrelevant since they 
are almost the same and the experimental resolution 
does not allow us to distinguish the two. The four 
periodic states [19] share the same geometric mean 
L(l/2) = 3V^/{2V2)Lj = l.SALj. 

The calculated sheet inductances L{f) are theoretical 
predictions based on the XY model in absence of ther- 
mal fluctuations and so have to be compared with low 
temperature data taken at high frequency where the cor- 
responding time scale is too short to allow for fluctuations. 
Figure [H] shows ratios of L~^{f) extracted from low tem- 
perature inverse magnetoinductance isotherms like those 
in figH] but at much higher frequency (w/27r = 16kHz) 
and at / = 0, 1/3, 1/2. The data in figlS] unambiguously 
demonstrate that at low temperature the modulation of 
the measured sheet inductance is weaker than it would 
be in the related XY model in absence of thermal fluctua- 
tions, and the selected ratios tend to unity. This behaviour 
can definitively be ascribed to magnetic screening effects. 

Conclusion. — The field modulation of the sheet in- 
ductance in Josephson-junction arrays on a dice lattice 
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Fig. 5: (colour on-line) Low temperature sheet inductance ra- 
tios I/~^(/i)/L^^(/2) at selected frustration states /i,/2 = 
0,1/3,1/2 extracted from high frequency {uj/2n = 16kHz) 
magnetoinductance isotherms. Horizontal (coloured) dashed 
lines show the values of the ratios calculated in the correspond- 
ing frustrated XY ground states (GS). The arrow indicates the 
temperature at which A = a (extracted from figE]) 



is shown to vanish while lowering the temperature. The 
response is field-independent below a reduced tempera- 
ture T = ksT/Ej « 10^"^, or cquivalently below an easily 
accessible temperature T « 3K. Note that this threshold 
temperature may change in samples with substantially dif- 
ferent geometrical parameters. The very low temperature 
regime is characterised by a penetration length shorter 
than the lattice parameter, and the magnetic energy asso- 
ciated with screening currents is greater than the Joseph- 
son coupling energy. The observed behaviour is explained 
in terms of increasing magnetic screening, with decreasing 
temperature, until a strong screening regime is reached 
where the modulation of the measured sheet inductance 
is weaker than in the frustrated XY model in absence of 
thermal fluctuations. The disappearance of the modula- 
tion in the magnetoinductance at low temperature could 
also be attributed to vortices with very low mobility. How- 
ever, this would give rise to a hysteretic response that was 
not observed down to the lowest accessible temperatures. 
In the (very) low temperature screening regime, the JJA 
behaves as a multiply connected superconductor, and the 
XY model is no longer valid. Unlike the weak screening 
regime where vortices are phase configurations satisfying 
the fiuxoid quantisation, in this strong screening regime 
vortices are real magnetic objects carrying an integer num- 
ber of flux quanta and interacting with currents through 
Lorentz forces. These strong effects highlight the impor- 
tance of a careful interpretation of low temperature results 
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within the frustrated XY model. 
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